Effects of turbulent mixing on critical behaviour: 
Renormalization group analysis of the Potts model 

N V Antonov and A V Malyshev 

Department of Theoretical Physics, St. Petersburg State University 
Uljanovskaja 1, St. Petersburg, Petrodvorez, 198504 Russia 

E-mail: nikolai.antonov@pobox.spbu.ru, alvlamal@gmail.com 

Abstract. Critical behaviour of a system, subjected to strongly anisotropic turbulent 
mixing, is studied by means of the field theoretic renormalization group. Specifically, 
relaxational stochastic dynamics of a non-conserved multicomponent order parameter 
of the Ashkin-Teller-Potts model, coupled to a random velocity field with prescribed 
statistics, is considered. The velocity is taken Gaussian, white in time, with correlation 
function of the form oc S(t — t')/|kx| 1 , where k_L is the component of the wave 
vector, perpendicular to the distinguished direction ("direction of the flow") — the d- 
dimensional generalization of the ensemble introduced by Avellaneda and Majda [1990 
Commun. Math. Phys. 131 381] within the context of passive scalar advection. This 
model can describe a rich class of physical situations. It is shown that, depending on 
the values of parameters that define self-interaction of the order parameter and the 
relation between the exponent £ and the space dimension d, the system exhibits various 
types of large-scale scaling behaviour, associated with different infrared attractive 
fixed points of the renormalization-group equations. In addition to known asymptotic 
regimes (critical dynamics of the Potts model and passively advected field without 
self-interaction), existence of a new, non-equilibrium and strongly anisotropic, type 
of critical behaviour (universality class) is established, and the corresponding critical 
dimensions are calculated to the leading order of the double expansion in £ and e = 6 — d 
(one-loop approximation). The scaling appears strongly anisotropic in the sense that 
the critical dimensions related to the directions parallel and perpendicular to the flow 
are essentially different. 
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1. Introduction 

Numerous systems of very different physical nature reveal interesting singular behaviour 
in the vicinity of their critical points. Their correlation functions exhibit self-similar 
(scaling) behaviour with universal critical dimensions: they depend only on few 
global characteristics of the system (like symmetry or space dimension). Consistent 
qualitative and quantitative description of the critical behaviour is provided by the field 
theoretic renormalization group (RG). In the RG approach, possible types of critical 
regimes (universality classes) are associated with infrared (IR) attractive fixed points of 
renormalizable field theoretic models. 

Most typical equilibrium phase transitions belong to the universality class of 
the O ra -symmetric ip 4 model of an n-component scalar order parameter ip. Universal 
characteristics of the critical behaviour depend only on n and the space dimension d 
and can be calculated within systematic perturbation schemes, in particular, in the form 
of expansion in e = 4 — d, the deviation of the space dimension from its upper critical 
value d = 4; see the monographs [H [2] and the literature cited therein. 

Another important example is provided by the Ashkin-Teller-Potts (ATP) 
class of models [3] [9]. In the continuous formulation, they are described by the 
effective Hamiltonian for the n-component order parameter with a trilinear interaction 
term, invariant under the hypertetrahedron symmetry group [5] [9]. Such models 
have numerous physical applications: magnetic materials and solids with nontrivial 
symmetry, Edwards- Anderson spin-glass models within the replica formalism [TO], and 
so on. In general, the ATP models describe systems which locally have n states, but the 
energy of any given configuration depends on whether pairs of neighboring sites are in 
the same state or not [3] . The case n = 2 corresponds to nematic-to-isotropic transitions 
in the liquid crystals [IT] , while the formal limits n = and n = — 1 correspond 
to the percolation problem and the random resistor network, respectively [i~2]-|14j. 
Recently, models with trilinear interaction have attracted a new amount of interest due 
to their interesting formal properties [15j [16] and applications to the dynamics of first- 
order phase transitions [17]. The application of the cubic model to the Yang-Lee edge 
singularity has long been known [18] . 

The problem of the nature of the phase transition in the ATP model has a long 
and rather entangled history; see e.g. [I] [E] and references therein. According to 
Landau's phenomenological theory, existence of a trilinear term excludes the possibility 
of the second-order transition. On the contrary, exact two-dimensional results, numerical 
simulations and RG analysis suggest that for small n, the phase transition in the ATP 
model is of the second order, while for n large enough (n > 3 in two dimensions pE] 
and n > 7/3 in the vicinity of the upper critical dimension d = 6 [7J [8]) the transition 
becomes a first-order one. In this paper, we accept the point of view that the existence 
of an IR attractive fixed point of the RG equations implies the existence of a self-similar 
(scaling) asymptotic regime and thus the existence of a kind of critical state. 

It is well known that dynamical critical behaviour (critical singularities of relaxation 
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and correlation times, various kinetic and transport coefficients etc.) appears much 
richer, less universal and is comparatively less understood. Different nature of the order 
parameter (conserved or non-conserved), inclusion of additional slow modes (densities 
of entropy or energy) and interaction with hydrodynamical degrees of freedom produce 
different types of critical dynamics for the same static model j2j [19]. 

The behaviour of a real system near its critical point is extremely sensitive 
to external disturbances, gravity, geometry of the experimental setup, presence of 
impurities and so on; see the monograph [20J for the general discussion and the 
references. "Ideal" equilibrium critical behaviour of an infinite system can be 
obscured by limited accuracy of measuring the temperature, finite-size effects, finite 
time of evolution (ageing) and so on. In the presence of a distinguished direction, 
scaling behaviour of such systems can become strongly anisotropic, with different 
critical dimensions corresponding to different spatial directions. What is more, some 
disturbances (randomly distributed impurities in magnets and turbulent mixing of fluid 
systems) can change the type of the phase transition (first-order to second order one, 
and vice versa) and produce completely new types of critical behaviour (universality 
classes) with rich and rather exotic properties. 

Investigation of the effects of various kinds of deterministic or chaotic flows (laminar 
shear flows, turbulent convection and so on) on the behaviour of the critical fluids (like 
liquid crystals or binary mixtures near their consolution points) has shown that the flow 
can destroy the usual critical behaviour: it can change to the mean-field behaviour or, 
under some conditions, to a more complex behaviour described by new non-equilibrium 
universality classes [2T] [32]. 

In this paper we study effects of turbulent mixing on the dynamical critical 
behaviour of the systems, described by the generalized ATP model, paying special 
attention to anisotropy of the flow. Bearing in mind application to liquid crystals or 
percolation in liquid media, we consider a purely relaxational stochastic dynamics of a 
non-conserved order parameter of the ATP model, coupled to a random velocity field 
with prescribed Gaussian statistics. 

Recently, the models involving passive scalar fields advected by such synthetic 
velocity ensembles attracted a great deal of attention because of the insight they offer 
into the origin of intermittency and anomalous scaling in the real fluid turbulence; see 
the review paper [33] and references therein. The RG approach to the problem of 
passive advection is reviewed in [31]. In spite of their relative simplicity, such models 
reproduce many of the anomalous features of genuine turbulent heat or mass transport 
observed in experiments. In the context of our study, it is especially important that 
they allow to easily model anisotropy of the flow, which in more realistic models would 
be introduced by the initial and/or boundary conditions. More specifically, we employ 
the (i-dimensional generalization of a strongly anisotropic ensemble introduced in [35] 
and further discussed in a number of papers, e.g. [36] [38], in connection with the 
passive scalar problem: the velocity field is oriented along a chosen direction n and its 
correlation function depends only on the coordinates perpendicular to n. 
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Although simplified, the model appears rather nontrivial and captures the main 
property of the problem: existence of a new, non-equilibrium and strongly anisotropic, 
universality class of scaling behaviour. 

The plan of the paper is as follows. In section [2] we present the detailed 
description of the model and its field theoretic formulation. In section [3] we analyze 
the ultraviolet (UV) divergences, relaying upon the power counting and additional 
symmetry considerations. We show that the model, after proper extension, appears 
multiplicatively renormalizable. Thus we can derive the RG equations and introduce 
the RG functions ((3 functions and anomalous dimensions 7) in the standard manner; 
see section HI 

In section [5] we show that, depending on the relation between the spatial dimension 
d and the exponent £ in the velocity correlator, the model reveals four different types of 
critical behaviour, associated with four fixed points of the corresponding RG equations. 
Three fixed points correspond to known regimes: Gaussian or free field theory, non- 
interacting scalar field passively advected by the flow (the ATP nonlinearity in the 
original dynamical equations appears irrelevant), and the original critical behaviour of 
the model without mixing. The most interesting fourth point corresponds to a new full- 
scale non-equilibrium universality class, in which both the nonlinearity and turbulent 
mixing are relevant. 

The corresponding critical dimensions can be calculated as double expansions in two 
parameters: £ and e — 4 — d. The scaling behaviour appears strongly anisotropic in the 
sense that the critical dimensions related to the directions parallel and perpendicular 
to the flow are essentially different. The practical calculation of the renormalization 
constants, RG functions, regions of stability and critical dimensions was performed in 
the leading order (one- loop approximation); some of the results, however, are exact 
(valid to all orders of the double e-£ expansion). These issues are discussed in section 
El while section [7] is reserved for conclusion. 

2. Description of the model and the field theoretic formulation 

Relaxational dynamics of a non-conserved n-component order parameter (p a {x) with 
x = {t, x} is described by a stochastic differential equation 



where d t = d/dt, A is the (constant) kinetic coefficient, and rj a (x) is a Gaussian random 
noise with zero mean and the pair correlation function 

(T) a (x)r}b{x')) = 5 ab D v (x - x'), 




(2.1) 




(2.2) 
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d being the dimension of the x space. Near the critical point, the Hamiltonian ~H(<p) of 
the ATP model is taken in the form [5j El [7] 

n(<p)= J dx j-^ a (x)dV a (x) + -|y? a (x)y? a (x) + 

+ |y J R a 6 C V 9 a(x)v5 b (x)v9 c (x)| , (2.3) 

where di = d/dxi is the spatial derivative, d 2 = didi is the Laplacian, r oc (T — T c ) 
measures deviation of the temperature (or its analog) from the critical value and 
go is the coupling constant. Summations over repeated indices are always implied 
(a,b,c = l,...,n and i = l,...,d); after the functional differentiation in (12.1 j) one 
has to replace <^(x) — )■ 

Following |9j, we consider the generalized case of certain symmetry group Q, which 
has the only irreducible invariant third-rank tensor R a b c 'i with no loss of generality it is 
assumed to be symmetric. In the original ATP model Q is the symmetry group of the 
hypertetrahedron in n dimensions. Then the tensor R a bc is conveniently expressed in 
terms of the set of n + 1 vectors e a which define its vertices [3 [6] : 

Rabc — 1 



Ct Ct Ct 
o p p 

J a^b °c i 



where the e" satisfy 

n+1 n+1 n 

E e " = °> E« = (™ + 1 )^> E^ef = (" +1)^-1- ( 2 - 4 ) 

a=l a=l a=l 

Using equations ( 12 Ah all the contractions with the tensor R a bc can be calculated. For 
example, the following contractions of two and three tensors have the form 

RabcRabe = Rl^ce, RaecRchbRbfa = R2Rehf, (2.5) 

where 

R 1 = (n+ l) 2 (n- 1), i? 2 = (n + l) 2 (n-2). (2.6) 

Coupling with the velocity field v = {vi(x)} is introduced by the replacement 

d t ^Wt = d t + v i d i , (2.7) 

where V* is the Lagrangian (Galilean covariant) derivative. For incompressible fluid, 
the velocity field v is transverse due to the continuity relation: diVi = 0. The velocity 
ensemble is defined as follows [35]. Let n be a unit constant vector that determines some 
distinguished direction ('direction of the flow'). Then any vector can be decomposed 
into the components perpendicular and parallel to the flow, for example, x = x+ +nii 
with xi • n = 0. The velocity field will be taken in the form 

v = nt)(t,x+), 

where f(t,x+) is a scalar function independent of x\\. Then the incompressibility 
condition is automatically satisfied: 

d^i = d\\v(t,x ± ) = 0. 
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For v(t,x±) we assume a Gaussian distribution with zero mean and the pair 
correlation function of the form: 

/dk 
^-pexp{ik-(x-x') = 

j—^ exp{ik ± • (x ± -x' ± )}5„(fc ± ), fc ± = |kj (2.8) 

with the scalar coefficient functions 

D v (k) = 2n8(k\{) D v (k ± ), D v (k ± ) = D fc^ 1 "*. (2.9) 

Here and below d is the dimension of the x space, D > is a constant amplitude factor 
and £ an arbitrary exponent. The IR regularization in ( 12 .8p is provided by the cutoff 
k± > m (by dimension, r oc m 2 ). Precise form of the IR regularization is inessential; 
sharp cutoff is the most convenient choice from the calculational viewpoints. The natural 
interval for the exponent is < £ < 2, when the so-called effective eddy diffusivity has 
a finite limit for m — > 0; it includes the most realistic Kolmogorov value £ = 4/3. 

In order to ensure multiplicative renormalizability of the model, it is necessary to 
split the Laplacian in (12. ip into the parallel and perpendicular parts d 2 — > dj_ + f d 2 
by introducing a new parameter fo > 0. Here d\ is the Laplacian in the subspace 
orthogonal to the vector n and dn = d/dxn. In the anisotropic case, these two terms 
will be renormalized in a different way; more detailed discussion of this point can be 
found in [3TJ [32] . Thus equation ( 12. ip becomes 

V t <^a = A (d 2 ± + f d 2 - T ) (f a ~ R a bc Vb^c + Va] (2.10) 

this completes formulation of the model. 

According to the general rule (see e.g. chap. 5 of the monograph [2J), our 
stochastic problem is equivalent to the field theoretic model of the extended set of 
fields $ = {ip' a , ip a , v} with the action functional 

Sm = + V'a {-V, + A (d 2 ± + fodl) - A ro} <p a - 

,1/4 

RabcV'aVbVc + S v (v), (2.11) 



Aofi'o/o 



where we segregated the factor / 1//4 from the charge go- The first few terms represent 
the De Dominicis-Janssen action functional for the stochastic problem (12. ip . (I2.2p at 
fixed v; it involves the auxiliary scalar response field <f' a {x). All the required integrations 
over x = {t, x} and summations over the vector indices are implied, for example, 

<Padi.<Pa = ^2 dt dx(p' a (x)d 2 ± (p a (x). 

0=1 J J 

The last term in (12.111) corresponds to the Gaussian averaging over v with the 
correlator (12. 8p and has the form 

S v {v) = ~ j ' dt J dx ± dx.' ± v(t,x ± )D^ 1 (yi ± -x ± )v(t,x A _), 
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9o 
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d F 
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1 

















dp 


1/2 


1/2 


-1 





-2 














d F 


(d + l)/2 


(d-3)/2 





-2 


2 


2 


1 


(6 - d)/2 





where 

is the kernel of the inverse linear operation D~ x for the correlation function D v in (12. 9ft . 

This formulation means that statistical averages of random quantities in the original 
stochastic problem coincide with the Green functions of the field theoretic model 
with the action (12. lip , given by functional averages with the weight exp<S(<&). This 
allows one to apply the standard Feynman diagrammatic technique, the field theoretic 
renormalization theory and renormalization group to our stochastic problem. 

3. Canonical dimensions, UV divergences and renormalization 

It is well known that the analysis of UV divergences is based on the analysis of 
canonical dimensions; see e.g. [TJ [2]. In general, dynamic models have two scales: 
canonical dimension of some quantity F (a field or a parameter in the action functional) 
is completely characterized by two numbers, the frequency dimension d F and the 
momentum dimension dp; see e.g. chap. 5 in [2J. They are determined such that 
[F] ~ [T] ~~ d F [L] ~ d F , where L is some length scale and T is the time scale. 

Our strongly anisotropic model, however, has two independent momentum scales, 
related to the directions perpendicular and parallel to the vector n, and requires a more 
detailed specification of the canonical dimensions. Namely, one has to introduce two 
independent momentum canonical dimensions dp and d F so that 

[F] ~ [T]-^[Lj-^[L||]- d K 

where L± and Lii are (independent) length scales in the corresponding subspaces. The 
dimensions are found from the obvious normalization conditions = —d^ ± = 1, 
d\ ± = —dx ± = 0, d% ± = d% = 0, d 1 ^ = —df = 1, and so on, and from the requirement 
that each term of the action functional (12.111) be dimensionless (with respect to all the 
three independent dimensions separately). 

The canonical dimensions of the model (12.111) are given in table [TJ including 
renormalized parameters, which will be introduced a bit later. From table Q] it follows 
that the model is logarithmic (the coupling constants go and Wq are simultaneously 
dimensionless) at d — 6 and £ = 0, so that the UV divergences in the correlation 
functions manifest themselves as poles in e = 6 — d, £ and their linear combinations. 

The total canonical dimension can be found from the relation dp = dp + d F +2d F (in 
the free theory, dt oc d\ oc <9jj). In the renormalization theory it plays the same part as 
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the conventional (momentum) dimension does in one-scale static problems: superficial 
UV divergences, whose removal requires counterterms, can be present only in those 1- 
irreducible Green functions, for which the total canonical dimension at the logarithmic 
values e = £ = (formal index of divergence) is a non-negative integer. 

The careful analysis of the table (TJ augmented by symmetry considerations, shows 
that all the counterterms needed to cancel the UV divergences in our model are 
present in the action (12.111) . Here, important role is played by the Galilean symmetry 
and the invariance with respect to the group Q. For example, the function ((p'(pvv) 
can be omitted from consideration because the corresponding counterterm (p' a (p a v 2 is 
forbidden by the Galilean invariance. A similar situation occurs for the functions 
((p'v) and (ip'vv). For the first function, a counterterm necessarily contains a spatial 
derivative, whose vector index is contracted with the velocity field; this gives zero 
in view of the transversality condition (12. 8ft . For the second function, the possible 
counterterms (p' a (diVk)(dkVi) and ip' a (diVk)(diVk) are forbidden by the symmetry with 
respect to Q. Furthermore, the Galilean symmetry requires that the counterterms 
ip' a dt^p a and <f' a (vidi)(p a enter the renormalized action only in the form of the Lagrangian 
derivative yj^VtW in turn, owing to the symmetry with respect to Q, the trilinear term 
in (12. lip is renormalized as a single entity. 

Thus our model appears multiplicatively renormalizable with the renormalized 
action of the form 

s R (a>) = Zl ^^ + A {-z 2 v t + a (z 3 dl + zjdf) - z 5 \r} cp a - 

\gn £ / 2 f 1 / 4 

~Z G Rabc<Pa ( Pb<Pc + Sv(v). (3.1) 

Here A, r, /, g and w are renormalized analogs of the bare parameters (with the 
subscripts "0") and fi is the reference mass scale (additional arbitrary parameter of 
the renormalized theory). The renormalization constants capture all the divergences at 
e, £ = so that the correlation functions of the renormalized model (13. ip have finite 
limits for e, £ — > 0, when expressed in renormalized parameters A, /i and so on. 

Simple analysis shows that the Feynman diagrams needed for the calculation of 
the renormalization constants in (13. ip are the same as in the one-component ip 3 theory 
(model 2 from [32]), multiplied with the appropriate tensor contractions. Moreover, in 
the one-loop approximation only contractions (12. 5 p appear. Therefore, one can easily 
generalize the results of [32] to the case at hand with the aid of (12. 5p . 

In [32] the minimal subtraction (MS) renormalization scheme was employed. In the 
MS scheme the renormalization constants have the forms Z$ = 1+ only singularities 
in e and £, with the coefficients depending on the two completely dimensionless 
parameters — renormalized coupling constants g and w. The one- loop results for the 
renormalization constants in (13. ip are the following 

12 2e ' d 3e ' 4 3e £ ' 

2uR 1 2uR 2 , , 

2-5 = 1 -, Z B = 1 -. (3.2) 
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Here we have passed to more convenient coupling constants u — > g 2 /128ir 3 and 

w -> w/24:7T 3 . 

The parameters R\ and R2 are related to the dimension n of the order parameter 
by the expression (I2.6p . Although we are especially interested in the cases n = and 2, 
for completeness the coefficients R\ and R2 in what follows are assumed to be arbitrary. 

Expression (I3.ip is equivalent to the multiplicative renormalization of the fields 
¥a -> PaZ^, <f r a —> (p'aZpi and the parameters: 

A = XZ X , r = r Z T , f Q = fZ f , 

g = gfi e/2 Z g (u = ufi £ Z u ) , w = wfJZ w (3.3) 

(no renormalization of the velocity field is needed: Z v — 1). The constants in Eqs. ( 13. ip 
and (13 .3p are related as follows: 

Z\ = Z\Z 2 ,, Z2 = Z^/Zp, Z3 = Z^/ZxZp, 

Z4 = Z^rZxZfZp, Z 5 = Z lp /ZxZ T Z lp , Zq = ZxZl/ 2 Z l J 4: Z ip iZ 2 p . (3.4) 

Since the last term <S„(v) is not renormalized, the amplitude D is expressed in 
renormalized parameters as 

A) = wo/oAo = wftfX, 
which leads to the relation 

Z w Z f Z x = 1 (3.5) 
between the renormalization constants. 



4. RG equations and RG functions 

Let us recall an elementary derivation of the RG equations; more detailed discussion 
can be found in monographs [fl [2]. The RG equations are written for the renormalized 
correlation functions Gr = (<&••• $)_r, which differ from the original (unrenormalized) 
ones G = ($•••$) only by normalization and choice of parameters, and therefore can 
equally be used for analyzing the critical behaviour. The relation Sr(Z$$, e, //) = 
5($,eo) between the functionals (12. lip and (13.11) results in the relations 

G{e ,...) = ZpZpG R {e,n,...) (4.1) 

between the correlation functions. Here, N v r and are the full numbers of 
corresponding fields entering into V (we recall that in our model Z v = 1); e$ = 
{Ao, To, fo, uo, wo} is the full set of bare parameters and e = {A, r, /, u, w} are 
their renormalized counterparts; the ellipsis stands for the other arguments (times, 
coordinates, momenta etc.). 

We use to denote the differential operation \xd^ for fixed e and operate on both 
sides of the equation (14. II) with it. This gives the basic RG differential equation: 

{D RG + N 9lv + iVv) G R (e, //,...) = 0, (4.2) 



Effects of turbulent mixing on the critical behaviour 10 

where T>rq is the operation D M expressed in the renormalized variables: 

V RG = P M + Mu + - IfDf - i x V x - 1t V t . (4.3) 

Here we have written T> x = xd x for any variable x, the anomalous dimensions 7 are 
defined as 

Jf = V^\nZ F for any quantity F, (4.4) 

and the j3 functions for the two dimensionless couplings u and w are 

j 14 = D Al u = u[-e-7 1t ], /Su, = V fl w = iu[-f-7«,], (4.5) 

where the second equalities come from the definitions and the relations (13. 3p . 

Equations H3.4f) result in the following relations between the anomalous dimensions 

7i = 7 A + 27^, 72 = 7<//+7</>, 73 = l v ' + 7a + l<p, 

74 = 7/ + 73, 75 = 7r + 73 , 76 = 7a + 7«/2 + 7//4 + 7^' + 27 v , (4.6) 
while from (I3.5P one obtains 

7/ + lw + 7a = 0. (4.7) 

The dimensions 71-76 are calculated from the corresponding renormalization constants 
using the definition (I4.4p : 

uR\ uRi uRi 

7i = 72 = -g-, 73 = 74 = -y + 

75 = 2uR u 7e = 2 Mj R 2 , (4.8) 

with the corrections of the order u 2 , w 2 , uw and higher. 

The RG functions entering the equation (14.31) are easily found from the relations 
(E2D, (TMD and (TCZD: 

uR\ 

Iff = 72 - 73/2 = uKi/3, 7^ = 73/2 = 

7/ = 74 - 73 = 7r = 75 - 73 = 5wi?i/3, 

7a = 73 - 72 = -uRt/6, 7«, = 72 - 74 = M-R1/6 - w, 

7, = 2 76 - 573/2 - 74/2 = (4i? 2 - u - w/2, (4.9) 

with the higher order corrections. 

From the definitions (I4.5j) . relations (14. 9p and explicit expressions H4.8I) for the 
anomalous dimensions we derive the following leading-order expressions for the (3 
functions of our model: 

(3 u = u [-e + Ru + w/2], W = w [-£ - uRj6 + w] , (4.10) 

where we have introduced a new convenient parameter R = Ri — 4R 2 . 
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5. Fixed points and scaling regimes 

It is well known that possible large-scale scaling regimes of a renormalizable model are 
associated with IR attractive fixed points of the corresponding RG equations. In our 
model, the coordinates w*, w* of the fixed points are found from the equations 

P u (u*, w*) = 0, f3 w (u*, u>*) = 0, (5.1) 

with the (3 functions given in (14. 5p . The type of a fixed point is determined by the 
matrix 

£1 = {% = dfii/du,}, (5.2) 

where A denotes the full set of the (3 functions and Uj = {u,w} is the full set of 
couplings. For IR stable fixed points the matrix Q is positive, i.e., the real parts of all 
its eigenvalues are positive. This condition defines the regions of IR stability for the 
corresponding scaling regimes. 

The couplings u and w should be non-negative (by definition, u cx g 2 > and 
w oc D /f\ > 0), and in the following we will be interested only in the "good" 
(admissible from the physics viewpoints) fixed points, which satisfy the conditions 

u* > 0, w* > (5.3) 

and can be IR attractive for some values of the model parameters. 

In order to give the complete picture of possible scaling regimes, it is instructive to 
discuss at first a more general case, specified by the (3 functions of the form 

A = u[— e + au + bw], (3 W — w[— £ + cu + dw], (5.4) 

with arbitrary real coefficients a-d. 

From equations (15. ip and (15. 4p we can identify four different fixed points. For 
the first three points the matrix Q appears to be triangular, so that its eigenvalues 
(and hence the regions of IR stability of the corresponding scaling regimes) are simply 
determined by the diagonal elements: Q u = d(3 u /du > and Q w = d(3 w /dw > 0. 

1. Gaussian (free) fixed point: it* = w* = 0; Q u = —s, Q w = — £. 

2. u* — (exact result to all orders), = £/d; Q u = —e + b^/d, Q w = £. This 
point can be "good" only if d > 0; otherwise the the conditions fl w > and u>* > 
cannot be simultaneously satisfied. 

3. w* = (exact result to all orders), = e/a; Q u = e, Q w = — £ + ce/a. Similarly 
to the case 2, this point can be "good" only if a > 0. 

The last, fully nontrivial, fixed point requires more detailed discussion. 
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4. The coordinates of this point are 

= (efe — 6£)/A, = (a£ — ce)/A, A = ad — be, (5.5) 
while the matrix Q can be written in the form 
' au* bu* ' 

ft = . (5.6) 

l cw* dw* J 

(it is useful not to substitute explicit expressions (15. 5p into (15. 6p for a while). 

The necessary and sufficient condition for the IR stability of this point can be 
restated as the fulfillment of two inequalities: 

detft>0, trft>0. (5.7) 

From ( 15. 6 p one obtains 

det ft = ivu>* A > 0, (5.8) 

which along with (15. 3p shows that this point can be "good" only if A > 0. 
For the trace of ft we obtain: 

trft = cm* + dw* > 0. (5-9) 

There are three possibilities: 

1) a > 0, d > 0. In this case the inequality ( 15.91) is an automatic consequence 
of (15. 3p . Four regions of stability of the fixed points 1-4 divide the e— £ plane without 
"gaps" or overlaps. This is the most typical situation, realized for the ip A model or the 
Gribov process in various kinds of random flows [30]- [32]. 

2) a < 0, d < 0. Then ( 15. 9ft contradicts to ( 15.31) and this point can never be "good." 
Such a situation was not yet encountered. 

3) The parameters a and d are opposite in sign: ad < 0. For definiteness, we 
assume that a < 0, d > 0. We will see in short that this situation can realize for the 
ATP model. In this case one obtains from (I5.9P : 

> —au*/d > 0, (5.10) 

where the last inequality follows from a < 0, > 0. The second inequality in ( 15.31) is 
implied by the (15.101) and thus becomes superfluous. The region where the fixed point 
is IR attractive and positive is given by the two inequalities 

> 0, au* + dw* > 0. (5-11) 

Summing up, we conclude that the scaling regime corresponding to the fixed point 
4 exists if A > and at least one of the two parameters a and d is positive. The region 
where the fourth fixed point is "good" is determined by the inequalities (15. 3p if a and 
d are simultaneously positive, and by the conditions of the type (15. 1 IR if a and d are 
opposite in sign. 
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Let us turn to our specific model with the (3 functions (14. 10p . Identifying them 
with ( 15. 4 p gives 



The coordinates of the four possible fixed points are obtained by substituting the 
expressions ( I5.12p into the general results. 

In the scaling regime corresponding to the fixed point 2, the nonlinearity if 2 in 
the stochastic equation (12. ip becomes irrelevant in the sense of Wilson due to the 
exact relation u* = 0. Thus we arrive at the linear advection-diffusion equation for 
a passive scalar field ip. In turn, the effects of the velocity field become irrelevant 
in the third regime (fixed point 3) due to the exact relation w* = 0. The isotropy 
violated by the velocity ensemble is restored and the leading terms of the IR behaviour 
coincide with those of the equilibrium dynamic model ATP. Finally, the fixed point 4 
corresponds to a new nontrivial IR scaling regime, in which the both nonlinearities in 
the stochastic equation for ip are important; the corresponding critical dimensions reveal 
strong anisotropy, depend essentially on the both RG expansion parameters e and £, 
and are calculated as double series in those parameters; see section El 

The regions of IR stability for all possible fixed points in the e-£ plane for different 
values of the parameters R\ and R2 are shown on figures 1-5. In the one-loop 
approximation, all the boundaries of the regions are given by straight lines. 

The figures [3] and [5] show that there are overlaps between the IR stability regions 
of fixed points 2 and 3 and points 1 and 4, respectively. This means that for the values 
of £ and £, corresponding to the region of an overlap, the system has two variants of the 
IR scaling behaviour. Which one of them is realized depends on the initial data for the 
parameters u and w in the RG equations. 

On the other hand, from figure H] one can see that if the parameters R\ and R2 
are such that Ar < and R < 0, there is a gap. The system does not exhibit scaling 
behaviour for corresponding values of e and £, which can be interpreted as existence of 
a first-order phase transition. 

It is interesting to note that for a < 0, the original static model has no "good" 
fixed point (which is usually interpreted as a first-order transition), but a "good" point 
of the type (I5.5P can appear in the full dynamic model with two couplings, as illustrated 
by figure 5. One can say that the phase transition changes its type and becomes a 
second-order one owing to the turbulent mixing. 

Let us briefly discuss the pattern of the RG flows in the plane of couplings u- 
w for two special cases. Consider first the simplest (and the most typical) situation, 
illustrated by figure [2j The schematic picture of the fixed points and RG flows is 
sketched on figure [6] for the case, when the values of the parameters e and £ lie in 
the region 4 on figure [2j Then all the fixed points lie in the physical domain u > 0, 
w > 0: the point 4 is IR attractive, point 1 is repulsive and points 2 and 3 are saddle 
points. If the parameters e and £ are changing such that the corresponding point in the 



a = R x - AR 2 = R, b = - 
A = R + R 1 /12 = A R . 



Ri 



d= 1; 




(5.12) 
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£ = 2e 




2 










4 










m 




S 


f = 


= -eRj / 6R 



Figure 1. A R > 0, R > 0, 
#1 > 0. 





Figure 2. A R > 0, i? > 0, 

i?i < 0. 




Figure 3. A R < 0, i? > 0, 
i?i < 0. 



Figure 4. A H < 0, R < 0, 
i?i arbitrary. 



f = 2e 




Figure 5. A fl > 0, i? < 0, i?i > 0. 
Regions of stability of the fixed points in the model ( 12. lip 
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W " 



W " 




Figure 6. RG flows for the 
case shown in figure 2. 



Figure 7. RG flows for the 
case shown in figure 3. 



e— £ plane crosses the boundary £ = 2e and moves into the region 2, the fixed point 4 
on figure 6 crosses the ray u = 0, w > 0, going through the point 2, and moves into 
the unphysical domain u < 0. The point 2 becomes IR attractive. If the point in 
the figure |2] crosses the boundary £ = —eRi/6R, moving from region 4 to 3, the fixed 
point 4 crosses the ray u > 0, to = 0, going through the point 3, and moves into the 
unphysical domain w < 0; the point 3 becomes IR attractive. Although this picture is 
based on the explicit one-loop expressions for the (3 functions and regions of stability, it 
appears robust with respect to the higher-order corrections. It is crucial here that the 
functions (3 U at w = and (3 W at u — coincide with the (3 functions of the corresponding 
single-charge problems: the ATP model and the passive scalar case, respectively. The 
latter are supposed to have a unique nontrivial fixed point each: for the passive scalar, 
the (5 function is known exactly, why for the ATP model (with R > 0) this is true at 
least within the e expansion. Thus we may conclude that the coincidence of the fixed 
points 2 and 4 or 3 and 4 takes place simultaneously with the changeover in their type of 
stability. In turn, this means that, although the boundaries between the regions 2 and 4 
or 3 and 4 in the £-£ plane can become curved beyond the one-loop approximation 
due to the higher-order corrections to the (3 functions in (15 .4p . no gaps nor overlaps 
can appear between those regions. This is equally true for the boundaries between the 
regions 1 and 2 or 1 and 3 in figure 2, but there the boundaries are not affected by the 
higher-order corrections due to the simple exact expressions for the eigenvalues Q u , w at 
the Gaussian fixed point. 

The RG flows for another interesting situation, illustrated by figure |3j are depicted 
on figure [7] for the case when the values of the parameters e and £ lie in the overlap of 
the regions 2 and 3. Then the both fixed points 2 and 3 are "good" and the asymptotic 
behaviour of a flow depends on the initial data for u and w. The point 4 is unphysical: 
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although it lies in the domain u > 0, w > 0, inspection of the explicit expressions (15. 5|) 
and (15 .8p shows that it is a saddle point for this case. It is also worth noting that, for 
all possible situations, the RG flow with the initial data in the physical domain u > 0, 
w > can never leave it, because the function (3 U vanishes for u = and arbitrary w, 
while (3 W vanishes for w = and arbitrary u. 

Let us conclude this section with a brief discussion of the most interesting 
physical cases of the original ATP model. Then from (I2.6P and (I5.12p we obtain 
R={n + 1) 2 (7 - 3n) and A R = (n + 1) 2 (83 - 35n)/12. 

The case n = 2 corresponds to the nematic-to- isotropic transition in a liquid crystal. 
Then 

^ = ^1 = 9, R 2 = 0, A R = 39/4, (5.13) 

and the case represented by figure Q] is realized. One can see that the interval of the 
most realistic values of the model parameters, d = 3 and < £ < 2, belongs completely 
to the region of stability of the most nontrivial fixed point 4. 

The second case of interest is the critical behaviour of bond percolation, obtained 
in the limit n = 0. Strictly speaking, more realistic description of dynamical percolation 
is given by a special model with a nonlocal in-time interaction [39], but it is interesting 
to look at the dynamics of the ATP model with n = as some kind of approximation. 
For n = relations (12.61) give 

R 1 = -l, R 2 = -2, R = 7, A R = 83/12. (5.14) 

Thus, we arrive at the case shown on figure EJ Again, the most realistic values of 
the model parameters (d = 3 and £ = 4/3) belong to the stability region of the new 
anisotropic scaling regime, corresponding to the fixed point 4. 

For n > 3 the case shown on figure H] is realized. Now the values d = 3 and 
£ = 4/3 lie in the "empty" region where none of the fixed points are "good." This fact 
is usually interpreted as a first-order transition; the turbulent mixing does not change 
its type. The interesting situation illustrated by figure [5J where the mixing gives rise to 
the changeover in the type of the phase transition to the second-order one, is realized 
for the interval 7/3 < n < 83/35, which does not contain integer values and has no 
documented physical interpretation. 



6. Critical scaling and critical dimensions 

We recall the definition of generalized homogeneity. Let F be a function of n independent 
arguments {xi, ...,x n } that satisfies the dimensional relation 

F(X^x u X a -x n ) = X aF F( Xl , ...,x n ) (6.1) 

with a certain set of constant coefficients (scaling dimensions) {cti, . . . , a n , ap} and an 
arbitrary positive parameter A > 0. Differentiating relation ( 16. II) with respect to A and 
then setting A = 1, we obtain a first-order differential equation with constant coefficients 

n 

a(Di F(xi, ...,x n ) = a F F(xi, ...,x n ), Vi = Xid/dxi. (6.2) 

8=1 
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Its general solution has the form 



F(xi,x 2 , ■ ..,x n ) = x x F 



012 /on ' • • • ' a n /ai 

where F is an arbitrary function of (n — 1) arguments. Obviously, the dimensions 
are defined up to a common factor (this can be seen by replacing A — > X a in (16. ip or 
multiplying equation (16. 2p by a); this arbitrariness can be eliminated, for example, if 
we set a.i = 1. If a, = for some x iy this variable is not dilated in (16. ip . and the 
corresponding derivative is absent from (16. 2p . 

It is well known that the leading terms, determining the asymptotic behaviour of 
(renormalized) correlation functions at large distances, satisfy the RG equation (14. 2p . 
in which the renormalized coupling constants are replaced with their values at the fixed 
points. In our case, this leads to the equation 

- 1 )d s - 7*d a - 7 *d t + 7 ;n v + 7 ;<M GR = °> ( 6 - 3 ) 

where 7* = 7(14 = it*, w = w*) for all the anomalous dimensions. 

We are interested in the critical scaling behaviour, that is, behaviour of the type 
(16. ip in which all the IR relevant parameters (momenta/coordinates, frequencies/times, 
deviation of the temperature from its critical value r oc (T—T c )) are dilated, while the IR 
irrelevant parameters (those which remain finite at the fixed point: A, /1 and /) are fixed 
[U [2] . Thus we combine the equation (I6.3P with the analogous equations, corresponding 
to the canonical scale invariance (see section [3]), so that the derivatives with respect 
to the IR irrelevant parameters are eliminated; this gives the desired equation which 
describes the critical scaling behaviour (for more details, see e.g. [311 132]): 

{D ± + A||D,| + A^D U + A T D T - N„A„} G N<S> = 0, 

where D± = k±d/dk± and D» = k\\d/dh\. Here, = 1 is the normalization condition, 
and the critical dimension of any IR-relevant parameter F is given by the general 
expression 

A F = 4 + A||4 + A w d£ + 7£, (6.4) 

with canonical dimensions from Table [1] and the relations 

A w = 2- 7 *, A|| = (2 + 7 ;)/2. (6.5) 

We are in a position to write the final one-loop results for the critical dimensions. 
Substituting (14 .9p into the general formulae (16 .4p and (16.51) gives 

6 11 2 3 

d 1 u*Ri d u*Ri w* 

A 05 / = - + lH H , A w = H H . 

v 2 3 4' v 2 6 4 

By inserting the explicit expressions for the fixed point coordinates and taking the 

equality d = 6 — e into account one obtains the leading-order expressions for the critical 

dimensions. The results for all scaling regimes are summarized in Table 2. 
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The expressions for the first and second fixed points are exact. Other dimensions 
have corrections, given by higher powers of e for the third fixed point and higher powers 
of e and £ for the fourth one. The critical dimensions for the models of a liquid crystal 
and bond percolation are derived from the general results by substituting the expressions 
(15. 131) and ( 15. 14p . respectively. 

Let us discuss the consequences of the general scaling relations for the most 
interesting special case of the pair correlation function. They result in the scaling 
expression 

(^(x + r.i + fj^f)) = 5 ab r~ 2A - t(t q t^, t/r*», r\\/rty , (6.6) 

where r_|_ = |r_J, rii = |r||| and J 7 is some scaling function. This representation is valid 
in the symmetric phase (ro > 0), where the tensor structure is simply given by the S 
symbol. It is natural to assume that T has a finite limit for tq = (that is, exactly at 
the critical point) and/or for t = (equal-time correlation function). Then from (16.61) 
one obtains 

(V9 a (x + r,t) ip b (x,t)) = S^rf 2 ^ T fry / r^" \ 

with another nontrivial function T{x) = .F(0, 0,x). 

The two last arguments in the scaling representation (I6.6P can also be chosen in 
the form rj_/Lj_(t) and r\\/L\\(t) with two different characteristic length scales 

L ± (t)~t a \ L||(i)~ru, a ± = l/A w , ay^Aii/A,. (6.7) 

For the most realistic values e — 3 (d — 3) and £ = 4/3 (Kolmogorov spectrum of the 
velocity) and for the n = 2 case of our model (liquid crystals) the explicit results from 
Table 2 and the expressions ( 15. 13ft give 

A w ps 2.359, Ay « 1.846, a± « 0.424, a\\ ps 0.783, A v ps 1.487, 

while for the percolation limit n = from ( 15.141) one obtains 

A w ps 1.944, Ay Pd 1.639, a x ps 0.514, a\\ ps 0.843, A v ps 0.731. 

Existence of two different length scales (16. 7p with power-law dependence on the time was 
established earlier in a number of studies within numerical simulations [26], approximate 
analytical solutions [27] , RG analysis [3TJ [32] and exactly soluble simplified models [28J . 
It is interesting to note that the inequality aqi > a± also holds for all those cases. 
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7. Conclusion 

We studied effects of turbulent mixing on the critical behaviour of the system, described 
by a relaxational dynamics of a non-conserved order parameter of the ATP model. 
The mixing was modelled by a Gaussian statistics with vanishing correlation time and 
strongly anisotropic correlation function oc 8(t — t')/k^ 1+ ^; see equations (12.8p . (12. 9ft . 
Such ensembles were employed earlier in [35]-|37j in the analysis of the two-dimensional 
passive turbulent advection (linear equation for the scalar field). 

The model, originally described by stochastic differential equations (I2.1l) - (l2.3p . 
(12. 7j) . can be reformulated as a multiplicatively renormalizable field theory with the 
action (12.111) . which allows one to employ the field theoretic RG to study its critical 
behaviour. The model reveals four different IR scaling regimes, related with the four 
different fixed points of the RG equations. Their regions of stability in the £-£ plane 
were identified in the leading order of the double expansion in e and £ and are shown 
on figures 1-5. These regimes correspond to: 

(1) Gaussian (free) model; 

(2) Linear passive scalar advection (the self-action term in the ATP Hamiltonian 
(I2.3p is irrelevant in the sense of Wilson) ; 

(3) Equilibrium critical dynamics of the ATP model (interaction with the velocity 
field is irrelevant); and 

(4) The full-fledged strongly anisotropic scaling regime in which the both 
interactions are important; it corresponds to a new non-equilibrium universality class. 

It was shown that the equilibrium critical regimes for the both physically interesting 
cases (liquid crystals and percolation process) become unstable for the realistic range of 
parameters d = 3 and < £ < 2, which includes the Kolmogorov spectrum (£ = 4/3) 
and the Batchelor limit (£ — > 2) and is replaced with the new non-equilibrium regime. 
The corresponding critical dimensions were calculated to first order of the corresponding 
RG expansion, which in this case takes on the form of the double expansion in e and £; 
explicit expressions are given in Table 2. 

Those results were derived within the leading (one-loop) approximation, that is, 
in the leading order of the double expansion in e and £, and their validity for finite 
physical values of these parameters can be called in question (especially because of large 
physical values of e = 6 — d). Careful analysis of this problem requires calculation 
of the higher-order corrections and applying some kind of summation procedure to 
the results obtained, as was done e.g. in [9] for the scalar static ip 3 model. Such 
analysis goes far beyond the scope of the present paper, and we hope to address it in 
the future. Nevertheless, the discussion of the RG flows, given in section El suggests that 
the pattern of fixed points (and thus of critical regimes), obtained within the one-loop 
approximation, appears robust with respect to higher-order corrections and is preserved 
for finite values of e and £. 

Thus we hope that our simplified model of a non-conserved order parameter and 
Gaussian velocity ensemble captures the most important features of the full-fledged 
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problem: emergence of a new non-equilibrium universality class with a new set of 
critical exponents, completely different from those of the classical ATP model; existence 
(for a strongly anisotropic velocity ensemble) of two different length scales (with a 
power law time dependence), and so on. Further investigation should take into account 
conservation of the order parameter, compressibility of the fluid, non-Gaussian character 
and finite correlation time of the velocity field, and so on. This work is partly in progress 
and partly remains for the future. 
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